Spectrally Dispersed iC-Band Interferometric Observations of 
Herbig Ae/Be Sources: Inner Disk Temperature Profiles 



J.A. Eisner 1 - 2 , E.I. Chiang 1 ' 3 , B.F. Lane 4 - 5 , & R.L. Akeson 6 
j ae@astro . berkeley . edu 

ABSTRACT 

We use spectrally dispersed near-IR interferometry data to constrain the tem- 
perature profiles of sub-AU-sized regions of 11 Herbig Ae/Be sources. We find 
that a single-temperature ring does not reproduce the data well. Rather, models 
incorporating radial temperature gradients are preferred. These gradients may 
arise in a dusty disk, or may reflect separate gas and dust components with dif- 
ferent temperatures and spatial distributions. Comparison of our models with 
broadband spectral energy distributions suggests the latter explanation. The 
data support the view that the near-IR emission of Herbig Ae/Be sources arises 
from hot circumstellar dust and gas in sub-AU-sized disk regions. Intriguingly, 
our derived temperature gradients appear systematically steeper for disks around 
higher mass stars. It is not clear, however, whether this reflects trends in relative 
dust /gas contributions or gradients within individual components. 
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INTRODUCTION 



Disks of dust and gas around young stars provide laboratories for studying the initial 
conditions of planet formation. A knowledge of disk temperature profiles at stellocentric 
radii R < 1 AU is crucial for understanding both terrestrial and giant planet formation. 
Planet formation is only possible beyond the dust sublimation radius (where temperatures 
are < 1500 K), since dust is an essential building block of terrestrial planets and giant 
planet cores. The disk temperature profile also dictates the location of the "snowline" (e.g., 
Lecar et al.ll2006l ). beyond which temperatures are low enough for water ice formation. The 
snowline tells us where icy material originates, which is important for understanding how 
water-rich, Earth-like planets come to be. Furthermore, because the surface density in solids 
is higher at and beyond the snowline, giant planet formation is facilitated there. 

Spectral energy distributions (SEDs) provide some informati on about disk temperature 
profiles (e.g., iHillenbrand et al.l Il992l ; IChiang fc Goldreichl 119971 ). but cannot provide un- 
ambiguous constraints. The addition of spatially resolved information breaks degeneracies 
inherent in S EP modeling, and enables mapping of the density and temper ature structure 



of disks (e.g., lEisner et al.ll2005l ; Ivan Boekel et al.ll2004l ; IWilner et al.l 120001 ) . For example 



near-IR i nterferometric observations have shown that simple geom etrically thin disk mod- 
els (e.g., iLynden-Bell fc Pringld Il974l ; lAdams. Lada. fc Shul 119871 ) generally do not fit the 



data obtained for protoplan etary disks, while the inclusion of puffed-up inner disk edges 



Dullemond et ajj 



Eisner et all 12004 120051 ) . 



20011 ) allows models to match the data well for most sources (e.g. 



To further constrain inner disk temperature profiles, spatially resolved observations 
at multiple wavelengths are necessary. Since dust (and gas) emission is temperature and 
wavelength dependent, different radial disk temperature profiles produce distinct spatial 
distributions of emission at different wavelengths. Thus, measurements of emission intensity 
and size as a function of wavelength can be used to determine inner disk temperature profiles. 

Spectrally dispersed interferometry observations also have the potential to probe the 
spatial distribution of spectral line emission relative to continuum emission. Most of the mass 
in protoplanetary disks is contained in gas, which emits primarily through molecular and 
atomic transitions. Investigating spectral line emission is therefore crucial for understanding 
disk dynamics and the mass assembly and migration of gas giant planets. 

Here we present observations that spatially resolve the .fT-band (Ao = 2.2 /zm; AA = 0.4 
fim) emission from protoplanetary disks around intermediate mass pre-main-sequence stars. 
These observations are also spectrally dispersed across 5 spectral channels (R = 25). We use 
these data to investigate the temperature profiles of disks at stellocentric radii less than 1 
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AU. We also explore whether disk models that include only emission from dust are sufficient 
to explain our observations or whether gaseous emission, including spectral line emission 
from hot water vapor or CO, needs to be considered. 



OBSERVATIONS AND CALIBRATION 



The Palomar Testbed Interferometer (PTI) is a long-baseli ne near-IR Michelson inter- 
ferometer located on Palomar Mountain near San Diego, CA (jColavita et al.lll999l ). PTI 
combines starlight pairwise from three 40-cm aperture telescopes using a Michelson beam 
combiner, and the resulting fringe visibilities provide a measure of the brightness distribution 
on the sky (via the van Cittert-Zernike theorem). PTI measures normali zed squared vi sibil- 
ities, V 2 , which provide unbiased estimates of the visibility amplitudes (1 Colavit aj 1 1 9 99l ) ; V 2 
is unity for point sources and smaller for resolved sources. 

We observed 14 Herbig Ae/Be (HAEBE) sources with PTI between May 2 002 and Jan 



uaiy 2004. Properties of the broadband K emission were reported previously in lEisner et al. 



(120031 . [2004J). Here we consider the subset of 11 objects that were spatially resolved (but 
not over-resolved) in the previous observations, and use the spectrally dispersed PTI data 
to examine trends in the visibilities with wavelength across a 5-channel "spectrum." The 
K-band (2.2 /j,m) measurements were dispersed into 5 spatially filtered spectral channels on 
an 85-m North- West (NW) baseline for all 11 objects, on an 86-m South- West (SW) baseline 
for 8, and on a 110-m North-South (NS) baseline for 5. The NW baseline is oriented 109° 
west of north and has a fringe spacing of ~ 5 mas, the SW baseline is 211° west of north 
with a ~ 5 mas fringe spacing, and the NS baseline is 160° west of north and has a fringe 
spacing of ~ 4 mas. 

Visibility data were measured in 130-second "scans," which each consist of 5 equal 
time blocks. V 2 for the scan is given by the mean of these blocks, and the measurement 
uncertainty is the standard deviation. Measured V 2 for the science targets were calibrated 
by observing sources of known angular size and parallax (see Eisner et al. 2004 for details of 
calibrator sources). These calibrators were unresolved by the interferometer, and were close 
(within ~ 10°) to the target sources on the sky; when calibrating a target source V 2 , data 
from different c alibrators were we ighted according to their proximity to the target in time 
and angle (e.g.. lBoden et al.lll998l ). 



Based on previous observations of binary objects, the uncertainty in calibrated, broad- 
band i^-continuum V 2 measurements for PTI is ~ 0.01 — 0.05. The data in the narrow-band 
spectral channels typically have similar uncertainties. We verified this by calculating the 
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standard deviation of all V 2 measurements in each spectral channel (calibrated as described 
above) for an unresolved "check star" whose K magnitude is similar to those of our targets. 
The standard deviations of the V 2 values for the check star in spectral channels centered at 
2.0, 2.1, 2.2, 2.3, and 2.4 fim are 0.05, 0.03, 0.02, 0.01, and 0.01, respectively. The larger 
values at the shortest wavelengths arise because absorption by atmospheric water vapor leads 
to lower fluxes. We adopt these values as the noise floors in each spectral channel for our 
science dataQ 

We also discard a small amount of data for which low photon counts prevent the mea- 
surement of any meaningful signal. For very low photo n counts the m easured V 2 , which are 



unbiased, may be > 1 or < because of large scatter (IColavital Il999l ). The discarded data 
are all from the shortest-wavelength channel, in which atmospheric water vapor absorption 
leads to lower photon counts. We eliminate all visibilities in the 2.0 /im channel from nights 
where the majority of the V 2 measurements in this channel are > 1 or < 0. 

We derive crude spectra from the observed count rates in our PTI data. We first 
normalize the photometry by dividing by the mean flux for all spectral channels, and then 
calibrate instrumental and atmospheric effects for this normalized photometry with the same 
stars used to calibrate the V 2 measurements. These main-sequence calibrator stars have 
spectral types ranging from F7 to AO, and nearly flat near-IR spectra (H — K < 0.2 based 
on 2MASS photometry). Since most calibrators are not colorless A0V stars, the calibrated 
fluxes for the target stars may have some residual slopes. However, this does not appear 
to be a significant effect. Figure [T] shows the fluxes for four of our target stars and one 
"check star," all calibrated with the same set of calibrators. The check star exhibits a 
nearly flat spectrum, as expected based on the near-IR colors (H — K ~ 0.1 from 2MASS 
photometry). Thus, our differential spectral calibration appears reliable to within ~ 10%. 
We estimate (conservatively) that relative (channel-to-channel) fluxes for our target objects 
have uncertainties of 20%. If we remove the continuum slopes of our targets, no significant 
flux variations are observed above the estimated 20% uncertainties. 

One of our sample objects, MW C 1080, has a nearby companion with 0''78 separation 



and AK = 2.70 (lEisner et al.ll2004j ). Since we do not know the color of the companion star 
across the K band, we do not attempt to subtract contributions of the companion from 
our flux and V 2 data. Given the flux ratio of the two stars this may introduce additional 
uncertainties of ~ 10% in the modeling for this source. None of our other targets have known 



1 The results discussed below do not change substantially if higher or lower noise floors are used. As 
expected, however, the uncertainties in fitted parameters are larger and the reduced \ 2 values of best-fit 
models are smaller if a higher noise floor is assumed (and vice versa for a lower noise floor). 
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comp anions which would influence the visibility data at the few percent level (jEisner et al 
2QM). 



3. MODELING AND RESULTS 

Before modeling the circumstellar excess emission from our sample, we account for the 
effects of flux from the unresolved central star on the measured V 2 in §3.11 In §3.2[ we fit a 
simple geometrical model to the V 2 data for each spectral channel individually; our intent 
is to provide a simple picture, with few assumptions, of how the emission size depends on 
wavelength. In §3.3} we then analyze as a whole our entire dataset of spectrally dispersed 
flux and V 2 measurements from 2.0-2.4 fim O, fitting disk models that account only for 
continuum emission to all of the data simultaneously. We consider a single-temperature ring 
model, a two-ring model, and a disk model for which the temperature decreases smoothly 
with radius. In §3.41 we investigate the possible effect of spectral line emission on the spec- 
trally dispersed data. Finally we compute broadband spectral energy distributions (SEDs) 
for our best-fit models and compare these to observed SEDs from the literature in §3.51 



3.1. Circumstellar-to-Stellar Flux Ratios 



In order to model the circumstellar emission, we must first estimate the contribution to 
the measu red fluxes and V 2 from th e unresolved central stars. Following previous investiga- 
tors (e.g.. iMillan-Gabet et al.ll200ll ). we use spectral decomposition to estimate the relative 
fluxes of the unresolved star and the circumstellar excess. We determine the stellar flux 
for our sample objects by fitting a blackbody to optical photometry (from the literature; 
see Eisner et al. 2004), and then extrapolating to the A-band. The excess flux is given 
by the difference between the observed flux and the estimated stellar flux. To convert the 
normalized 2.0-2.4 /zm fluxes in each channel (jQ into absolute fl ux units, we use continuum 
A-band photometric flux measurements from lEisner et al.l (120041 ) . Derived circumstellar-to- 
stellar flux ratios for each spectral channel are listed in Table [D Because the stellar flux 
declines precipitously with wavelength while the excess flux increases with wavelength, the 
circumstellar-to-stellar flux ratio increases steeply across the A-band. 



As discussed in Eisner et al.l ( 120041 ). uncertainties in the A-band photometric measure- 
ments and uncertainties in optical photometry lead to errors in the estimated circumstellar- 
to-stellar flux ratios of ~ 10-20%. In addition, uncertainties in the channel-to-channel flux 
measurements (£J21) introduce errors in the circumstellar-to-stellar flux ratios of ~ 10-20%. 
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While many of our sample objects are variable at near-IR wavelengths (e.g., ISkrutskie et al. 



19961 ). the photometric K-band measurements were nearly contemporaneous with our PTI 



observations, and thus source variability should not be a major source of uncertainty. Thus 
the total uncertainties in estimated circumstellar-to-stellar flux ratios are ~ 20-30%. These 
uncertainties lead to small errors in the fitted angular sizes because the near-IR emission 
from our sample is dominated by t he circumstellar co mponent; typical uncertainties on fitted 
parameters are a few percent (see Eisner et al.l 120041 ) . 



3.2. Emission Radius versus Wavelength 



Before considering the complete dataset, which consists of V 2 and fluxes measured in 
5 spectral channels from 2.0 to 2.4 /zm, we provide an estimate of angular size as measured 
in each spectral channel. We employ a simple model for the circumstellar excess emission, 
and fit this model to the V 2 measured in each spectral channel. Each source was observed 
multiple times over multiple baselines, and we fit all of the V 2 measurements in a given 
channel simultaneously. 

Our chosen model is a uniform ring of fixed fractional width. The normalized V 2 of a 
ring, viewed face-on, is given by 
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where 8- m /2 is the inner angular radius of the ring, 1 + / is the ratio of the outer ring radius 
to the inner ring radius, r uv = [u 2 + v 2 ) 1 ^ 2 is the "uv radius," u, v are the projected baseline 
lengths on the sky, and J% is a first-order Bessel function. We take / = 0.2 for all objects. 
The only free parameter in the model is 6- in . 

Inclined geometries are included in this model, and in those discu ssed below, using a 



trans f ormation to a circ ularly symmetric coordinate system (see e.g., iMillan-Gabet et al. 
200ll ; lEisner et al.l 120031 ) . We assume values for the inclination and po sition angle based 
on previous modeling of the broadband K emission ( lEisner et al.l 120041 ) . We verified that 
these values were reasonable for the spectrally dispersed data by fitting ring models where 
inclination and position angle (in addition to #i n ) were free parameters. For all objects, the 
best-fit position angles and inclinations were consistent within the la uncertainties across the 
5 spectral channels. Furthermore, t hese values were co mpatible with those derived previously 
from wide-band V 2 measurements (lEisner et al.ll2004l ). We therefore fix the inclinations and 
position angles in the modeling presented here. 



Fitted inner ring radii as functions of wavelength are shown in Figure [2j We also indicate 
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the corresponding linear radii computed using the distances assumed in Eisner et al.l ( 12004 ). 
Clearly, size increases with wavelength for most objects. Such behavior is expected for 
irradiated disks and/or accretion disks, where the temperature decreases with increasing 
stellocentric radius. It is not consistent with a single-temperature ring at fixed radius, since 
the angular size of such a ring is independent of wavelength. For several objects, the inferred 
sizes are not obviously increasing with wavelength. Rather, the size of the circumstellar 
emission is independent of A within uncertainties (e.g., V1295 Aql), or is not a monotonic 
function (e.g., V1685 Cyg). 



3.3. Circumstellar Emission Models 

To quantitatively investigate the radial temperature structure suggested by Figure [21 in 
this section we fit simple models to our V 2 and photometry data for all spectral channels 
simultaneously. Here we assume that the models emit only continuum emission; in §3.41 we 
investigate whether better fits are obtained when spectral line emission is included. For all 
models we fit only our V 2 and flux data measured from 2.0-2.4 /im; in §3.51 we compute the 
broadband SEDs predicted by our best-fit models and compare to observations. 



3.3.1. Single-Temperature Ring 

We first consider the uniform ring model discussed above, but now fit the model to the 
complete dataset all at once (instead of fitting each channel individually). In addition, we 
now consider the measured fluxes as well as V 2 . This necessitates the introduction of a ring 
temperature, T r ; ng , in order to fit the photometry data. If we assume for simplicity that the 
ring emits as a blackbody, the flux emitted by the ring is 

FAR) = 2 -^Pb u (T rmg ) Rl (f + f/2) « 2 -^p Bu (T rmg ) R 2 n f. (2) 

Here B u is the Pl anck function, d is the distance, and i is the inclination of the ring (d and i 



are assigned from Eisner et al.l 120041 ) . We determine the values of R in and T ring that provide 



the best fit to all of the measured V 2 and fluxes for a given source (Table [2]). 

It is not surprising that the reduced x 2 values of our best-fit ring models are > 1 (Table 
|2]), since the ring model produces the same size in each channel, contrary to what is observed 
in the data (see §3.2p . We illustrate the poorness of fit with an example in Figure [3l which 
shows the measured V 2 along with the prediction of the best-fit single-temperature ring 
model for AB Aur. While \ 2 is not far from unity for several objects, this is probably because 
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of large measurement errors and does not necessarily imply that the single-temperature ring 
model provides a good physical description of the data. 



3.3.2. Two Rings 

The dependence of source angular size on wavelength suggests that we should include 
temperature-dependent radial gradients in the models. While a single-temperature ring 
model possesses no such gradients, the addition of a second ring at different stellocentric 
radius and temperature will lead to observed source structure that is a function of wave- 
length. Specifically, the cooler, outer ring will contribute relatively more emission at longer 
wavelengths, and thus the angular size of the two-ring model will tend to increase with 
wavelength. 

We assume the rings have fractional widths of 0.2, and we fix the radius of the inner ring 
at 0.1 AU. So long as the inner ring is essentially unresolved by our observations, the exact 
choice of inner ring radius will not substantially alter the quality of the fits. We verified this 
by performing fits where the inner ring radius was left as a free parameter; reduced \ 2 values 
did not improve with this extra degree of freedom. Smaller assumed values of the inner 
ring radius will, however, yield higher inner ring temperatures. Conversely, larger values of 
the inner ring radius give lower inner ring temperatures. However substantially larger inner 
ring radii can be ruled out since they would be spatially resolved, and would thus lower the 
quality of the fits. 

The free parameters of the model are the temperature of the inner ring, and the radius 
and temperature of the outer ring. This is obviously a simple description of a potentially 
complex system, but will demonstrate whether a multi-component continuum model is suit- 
able for explaining our data. As above, the data fitted include V 2 and fluxes measured from 
2.0-2.4 fjm Qj2D- 

V 2 values predicted by best-fit two-ring models are shown in Figure HI Best-fit param- 
eters and reduced x 2 values are listed in Table [31 For objects that display a monotonic 
increase in size with wavelength (see Figure [2]) the reduced x 2 values for this model are 
lower than for the single ring model ( §3.3.11) . In contrast, for T Ori, V1295 Aql, V1685 Cyg, 
and AS 442, sources for which the size does not appear to increase with wavelength, the 
single-temperature ring model provides a lower reduced \ 2 ■ The high temperature of the 
inner ring for MWC 1080 is due to our assumed inner ring radius of 0.1 AU. A larger inner 
ring radius — which would lead to a lower fitted inner ring tempera ture — is compatible with 



our data; at the source distance (~ 1 kpc; see Eisner et al.ll2003l ). an inner ring radius as 
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large as ~ 0.5 AU would still produce unresolved emission as assumed in our model. 



3.3.3. Disk with Smooth Radial Temperature Gradient 

As a generalization of the two-ring model ( §3.3.2j> . we construct a disk model with a 
continuous temperature gradient. If the temperature decreases over a range of disk radii 
then the cooler material at larger radii may contribute relatively more emission at longer 
wavelengths. We consider a disk model with a power-law temperature profile, 

Here, R is the stellocentric radius, i? in is the inner radius of the disk, and T in is the temper- 
ature at R = Ri n . The model is described by three free parameters: i?i n , T in , and a. 

We assume the model extends from R [n to an outer radius Rout = 10 AU. The exact 
choice of R ont is not crucial since most of the 2 fim emission arises from sub-AU radii. 
However if R mt is 1 AU then the model will not produce substantially different angular 
sizes at different radii; rather the emission would appear to arise from a narrow ring at all 
wavelengths and the model would resemble the single-temperature ring model described in 
§3.3.11 To explain the dependence of angular size on wavelength, we require R out > 1 AU. 

We divide the disk up into annuli and compute the flux and visibilities for each annulus. 
We assume each annulus emits as a blackbody so that the flux in an annulus of infinitesimal 
width dR (observed at the Earth) is 

27T COS 1 

dF v [R) = ~^B V (T) R dR, (4) 

where B v is the Planck function, d is the distance, and i is the disk inclination. The normal- 
ized visibilities for an annulus extending from Ri to R 2 = Ri + dR are given by the difference 
of visibilities for uniform disks having radii equal to Ri and R 2 : 



V(R) X<1 



(R\ - R\) 



7rr„ 



2%r uv R 2 \ ( 2nr uv R 1 

ti2J\ — ri — — ti\J] 



Xd \ Xd 



(5) 



The sum of dF v over all annuli gives the total disk flux, F v . The square of the flux- weighted 
sum of V(R) over all annuli gives V 2 . 

The V 2 values predicted by best-fit models including a temperature gradient in the 
inner disk are shown in Figure HI Table H] shows that the reduced x 2 values of the best fits 
are generally < 1. For all objects except V1295 Aql and V1685 Cyg this model provides a 
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better fit (i.e., a lower reduced x 2 value) than the single-temperature ring model. Moreover, 
the temperature gradient model yields a better fit to the data than the two-ring model 
for all sources. Thus, our spectrally dispersed interferometry data favor disk models that 
incorporate radial temperature gradients. How well the temperature gradient model fits 
broadband SEDs compared to other models is discussed in §3.51 



For most so urces the inner dis k sizes and temperatures derived here are similar to 
those inferred by lEisner et al.l (120041 ) from broadband (2.0 -2.4 iim) V 2 measu rements. This 
is expected since the broadband V 2 measurements from lEisner et al.l (120041 ) have similar 
values to the spectrally dispersed measurements presented here. However, for objects with 
small best-fit a values (e.g., MWC 758 or CQ Tau), our determined inner radi i are smaller 



than t he previous estimates by < 40%. This discrepancy is due to the fact that lEisner et al. 
(120041 ) assumed a = 0.75 for their disk models; for the smaller a values inferred here 



the inner edge of the disk must extend further inward to match the data. The values of 
i?i n determined here are more reliable than previous estimates since a and R m have been 
determined simultaneously. 



3.4. Spectral Line Emission 

While disk models including only continuum emission can provide reasonable fits to our 
spectrally dispersed interferometry data, several objects — T Ori, V1295 Aql, V1685 Cyg, 
and AS 442 — do not show the monotonic increase of size with wavelength expected for the 
continuum models with temperature gradients (§ §3.3.2H3.3.3p . This is particularly striking 
for V1685 Cyg, where the smaller uncertainties strongly imply non-monotonic dependence of 
size on wavelength (see Figured]). One potential explanation is that gaseous emission, with 
a different spatial scale than the dust continuum, contributes to the visibilities in certain 
spectral channels. 

Based on compil ations of H 2 and CO emission lines from the HITRAN database 



( IRothman et al.l 120051 ) . gaseous spectral line emission, if present, contributes principally to 
the PTI measurements in the channels centered at 2.0, 2.3, and 2.4 /jm (Figure [5]). Previous 
observations of hot gas near a young star illustrate that gaseous emission features may be 



quite broad (jCarr et al.ll2004l ). and thus this emission might fill large parts of some spectral 
channels. 

If the dust and gas emission components have different spatial distributions, then the 
visibilities in spectral channels where gaseous emission is important (2.0. 2.3, and 2.4 /xm) 
may be offset from the visibilities in spectral channels dominated by dust continuum emission 
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(2.1 and 2.2 /im). This effect has been observed previously in evolved late-type stars, where 
molecular l ayers above the stellar photosphere lead to different measured sizes across the K- 



band (e.g.. iThompson et al.ll2002l ). To investigate whether gaseous emission is contributing 
to our measured visibilities, we fit the models discussed above to only the visibilities measured 
from 2.1 to 2.2 /im, and compare these fits to those using the data from all channels. 

Best-fit parameters for "continuum" -only fits (i.e., fits using only the data measured 
from 2.1 to 2.2 /im) are included in Table HI and the predicted V 2 values are plotted in 
Figure HI For all sources, the best-fit inner disk sizes for "continuum" channels only and 
for all channels are equal within the la uncertainties. While the reduced \ 2 values of the 
"continuum" -only fits are smaller for AB Aur, T Ori, VV Ser, AS 442, and MWC 1080 
(Table HJ), the differences are small for all sources except AB Aur. 

The incompatibility of a simple disk model with the observed size versus wavelength for 
VI 685 Cyg was our main impetus for investigating gaseous emission. However, even if we 
consider only channels that are not expected to be contaminated by gaseous emission, the 
data are still not fitted well by our disk models (note the large x? values in Table HJ); the 
circumstellar emission for this object may not be represented well by a pure disk model. 

The spectrally dispersed interferometry data presented here do not argue strongly for 
the presence of gaseous spectral line emission in the inner disks around our sample HAEBE 
stars. We also do not observe any strong gaseous emission features in our crude spectra (e.g., 
Figure [T]) . However the large uncertainties do not rule out gaseous spectral line emission 
with line-to-continuum ratios < 20%. Furthermore, the smaller Xr values of continuum-only 
fits for several objects (especially AB Aur) are intriguing. 



3.5. Spectral Energy Distributions 



We fitted the models described above to measured V 2 and fluxes measured from 2.0-2.4 
/im. We now compare the broadband fluxes predicted by our best-fit models to observed 
SEDs. T his comparison is motivated by previous ana lyses of SEDs and -fT-band interferome- 
try (e.g.. lMillan-Gabet et al.ll200ll ; Eisner et al.ll2004j ) . which showed that single-temperature 
ring models generally provided better fits to the data than disks with temperature gradients. 

Figure [7J shows the predicted and observed SEDs for our sample. Disk models incor- 
porating smooth temperature gradients do not fit the SED data as well as the two other 
models for almost all objects; the SEDs are fitted better by single-temperature ring models 
or by two-ring models. One exception is MWC 1080, for which a disk model fits the SED 
data well; for this source the two-ring model does not fit the SED well (Figure [7]). This is 
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consistent with previous arguments that early-type stars l ike MWC 1080 are surrounded by 



geometrically thin disks with temperature gradients (e.g.. Eisner et al.ll2004j ) 



4. DISCUSSION 

The size of the near-IR emission for most of our objects varies as a function of wave- 
length, indicating that a simple, single-temperature ring model for the emission is untenable; 
such a model predicts the same size regardless of observing wavelength. We showed above 
that the V 2 and flux data measured from 2.0-2.4 /im are fitted better by models incorporat- 
ing radial temperature gradients. 



Early puffed-up inner disk models (e.g.. iDullemond et al.ll200ll ) generally assumed that 



the inner wall emits as a single-temperature blackbody; this model resembled a single- 
temp erature ring. Even if temperature gradients within the puffed-up inner rim are included 
(e.g., Isella fc Nattall2005l ). the narrow range of radii over which the emission arises means 



that these models still resemble single-temperature rings. Our data therefore suggest that 
puffed-up inner rims alone do not fully represent the circumstellar emission from HAEBE 
sources. 

Temperature gradients over a range of disk radii must be included in the model to fit 
the observational data presented here. The best fits to our spectrally-dispersed interferom- 
etry are obtained for two-ring models or models with power-law disk temperature profiles. 
Comparison of our models with broadband SEDs suggests that two-ring models are prefer- 
able. We suggest that our two-ring model represents a puffed-up edge of the dust disk and 
hotter inner disk gas. The best-fit temperatures of the inner ring are substantially hotter 
than the dust sublimation temperature, supporting our hypothesis that such emission prob- 
ably traces gas. The outer rings have temperatures less than or comparable to the dust 
sublimation temperature, suggesting that the outer ring emission can arise from dust. 

If inner disk/ring temperatures larger than dust sublimation temperatures reflect con- 
tributions from hot gas emission to measured fluxes and V 2 , we might also expect to detect 
variations in the visibilities as a function of wavelength due to gaseous CO or H2O emission 
features. In general, no such trends are seen in our data ( §3.41) . It is possible that the gaseous 
emission features are so broad that they produce a "pseudo-continuum" ; there is some sup- 
port for this suggestion from p revious observations of broad gaseous emission features from 
a young star (ICarr et al.l 120041 ) . It is interesting to note, however, that for AB Aur, AS 442, 
and MWC 1080, objects with higher inner disk temperatures, there is a suggestion that the 
inclusion of gaseous emission features in the models would improve the fits to the observa- 
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tions, since the "continuum" -only models yield lower reduced \ 2 values than the fits where 
all the data is used ( §3.4p . Future observations with better spectral resolution are needed to 
probe potential spectral line emission in more detail. 

While a two-ring model provides the best match to the combination of spectrally dis- 
persed and SED data, a disk model with a radial temperature gradient provides the best 
fit to the spectrally dispersed interferometry data (from 2.0-2.4 /^m) for all sources. Our 
results suggest that a multi-component model that includes temperature gradients in one 
or more components may provide better fits to the data than the simple models considered 
here. We note that MWC 1080, the most massive star in our sample, provides an exception: 
disk models with smooth temperature gradients fit all of the data well. 

We discern a striking correlation between temperature gradient and effective tempera- 
ture of the host star T*, which may reflect relative differences in the dust/gas components or 
variations in the gradients within individual components. Using our power-law disk model, 
T = Ti n (R/Ri n )~ a , to illustrate this trend, we plot a as a function of T* in Figure El Larger 
values of a (steeper gradients) are associated with objects with earlier spectral types. For 
MWC 480, MWC 758, CQ Tau, and T Ori, the objects in our sample with the latest spectral 
types, the best-fit models require a = 0.40-0.55. In contrast, inferred values of a for other 
objects in our sample range between 0.60 and 1.0. 

There is some degeneracy between temperature profiles and optical depth profiles. In 
our modeling we assumed the emission was optically thick. The fitted inner disk and ring 
temperatures may be underestimated to the extent that the optical depth of the emitting 
region, t„, is < 1. By the same token, radial gradients in t v are somewhat degenerate with 
temperature gradients. However, since r v is likely to decrease with increasing radius and 
wavelength (e.g., r v oc R~ x for the surface layers of geometrically thin disks; Chiang & 
Goldreich 1997), r v gradients are more likely to produce larger sizes at shorter wavelengths 
(since more material at larger radii becomes optically thick). The gradients inferred from our 
data are therefore likely tracing temperature profiles as opposed to optical depth profiles. 

5. CONCLUSIONS 

We presented near-IR interferometry data, spectrally dispersed across 5 channels from 
2.0 to 2.4 fim, from the Palomar Testbed Interferometer. These data were used to constrain 
the temperature profiles of sub-AU-sized regions of 11 Herbig Ae/Be disks. We found that a 
single-temperature ring of emission did not fit the spectrally dispersed data well. Two rings 
at different stellocentric radii and with different temperatures provide better fits to the data. 
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The best fits (i.e., yielding the lowest reduced x 2 values) are for disk models incorporating 
smooth radial temperature gradients. 

However, the smooth gradient disk models do not match broadband SEDs nearly as 
well as two-ring models. We therefore speculate that multiple emission components, probably 
dust and gas, at different temperatures and stellocentric radii are contributing to the observed 
data. More complex two-component models, including temperature gradients in one or both 
components, will better match the observations. 

The inferred temperature profiles are systematically shallower for disks around less mas- 
sive stars. However, it is unclear whether this indicates relative differences in the properties 
of the two components (gas and dust) as a function of spectral type, or whether this trend 
traces temperature gradients within a single component. 

We also investigated whether gaseous emission, potentially at different spatial scales 
than the dust emission, contributes spectral line flux to our measurements. Because CO 
and H 2 are expected to emit primarily in the edge channels of our spectrally dispersed 
measurements, we performed fits of our disk models to only the central channels. Some 
objects showed a slight reduction (and in the case of AB Aur, a substantial reduction) 
in the reduced \ 2 values of the best-fit models when only the central channels were used. 
However, the best-fit parameters for fits to the central channels and those for fits to all 
data are consistent within the la uncertainties for all objects. We are therefore unable to 
claim strong evidence for spectral line emission in the inner disks of our sample based on our 
current data. Given the evidence for hot inner disk continuum emission, future observations 
with better spectral resolution and sensitivity are warranted to further investigate inner disk 
gas in these objects. 
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Table 1. Circumstellar-to-Stellar Flux Ratios from 2.0 — 2.4 /zm 



Source 


r 2.0/^m 


r 2.1 /Ltm 


r 2.2/^m 


r 2.3 fim 


r 2Afj,m 


AB Aur 


3.8 


4.7 


5.8 


7.0 


8.2 


MWC 480 


2.0 


2.6 


3.2 


3.9 


4.6 


MWC 758 


3.3 


4.1 


5.0 


6.0 


7.0 


CQ Tau 


2.5 


3.3 


4.2 


5.1 


6.1 


T Ori 


5.3 


6.7 


8.5 


10.3 


12.0 


MWC 120 


4.9 


6.3 


7.8 


9.5 


11.2 


VV Ser 


6.5 


8.2 


10.6 


13.1 


15.6 


V1295 Aql 


1.5 


2.0 


2.6 


3.3 


4.0 


V1685 Cyg 


3.7 


4.7 


6.2 


7.6 


9.1 


AS 442 


3.6 


4.7 


6.2 


7.7 


9.1 


MWC 1080 


5.2 


7.0 


9.1 


11.6 


14.1 



Note. — r\ denotes the ratio of the circumstellar excess 
emission to the stellar photospheric emission. This ratio 
is calculated for each of the PTI spectral channels, and 
typical uncertainties are ~ 20-30% ( §3.ip . 
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Table 2. Best fits for single-temperature ring models 



Source Xr (mas) R in (AU) T ring (K) 



AB Aur 


2.346 


i.64: 


hO.Ol 
-0.01 


23+ - 01 

u - zo -0.01 


1860.01 


MWC 480 


1.153 


1.511 


hO.Ol 
-0.01 


21 +o t)1 

U - Zi -0.01 


1440.0+ 


MWC 758 


2.267 


i.3o: 


hO.Ol 
-0.01 


20+ 01 
u -^ u -o.oi 


1600.01 


CQ Tau 


1.179 


i.36: 


h0.02 
-0.01 


20+ - 01 
u -^ u -o.oi 


1420.01 


T Ori 


0.904 


0.87: 


h0.02 
-0.03 


39+ - 01 

u - oy -0.01 


1640.01 


MWC 120 


1.173 


i.52: 


hO.Ol 
-0.01 


76+ - 01 

u - ' u -0.01 


1560.01 


VV Ser 


1.164 


i.44: 


h0.02 
-0.03 


45+ - 01 
u -^ d -o.oi 


1500.01 


V1295 Aql 


1.005 


L - ' u -0.01 


n 40+0.01 
u -^ y -o.oi 


1320.01 


V1685 Cyg 


4.101 


i.i6: 


1-0.01 
-0.01 


i - iU -o.oi 


1760.01 


AS 442 


1.125 


0.87: 


h0.03 
-0.04 


u - dz -0.02 


1840.01 


MWC 1080 


1.840 


i.34: 


hO.Ol 
-0.01 


i Q4+0.01 

i -°^-0.01 


2300.01 



-60.0 
-60.0 
-40.0 
-40.0 
-40.0 
-50.0 
-30.0 
-40.0 
-50.0 
-60.0 
-40.0 
-50.0 
-40.0 
-50.0 
-20.0 
-40.0 
-50.0 
-60.0 
-80.0 
-80.0 
90.0 
110.0 



Note. — Linear inner ring radii, R[ n , are computed using the 
de rived angular inner radii, 9- m /2, and the distances employed 
bv lEisner et aD (J20o3). 
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Table 3. Best fits for two-ring models 



Source Xr T im g ,i (K) 6» in>2 /2 (mas) R m , 2 (AU) T ring ,2 (K) 



AB Aur 


1.769 


2310.0™ 


o 1 4 +o.09 


S0+ 01 

u - ou -0.01 


1460.0±fg'g 


MWC 480 


0.841 


1710.0H85 


o i q+0.11 
z - ±o -0.16 


S0+ 02 

u - ou -0.02 


1130.0±lg|g 


MWC 758 


1.256 


2010.0lgg.-g 


2 24+0.01 

z - z ^-0.13 


S4+ - 01 

u -°^-0.02 


1060.0±|g-g 


CQ Tau 


1.107 


I 760 +80 

I I uu.u_ 140 


o on+0.04 

z - zu -0.41 


0.3333 


1010.0™ 


T Ori 


1.139 


31 00 n+ 690 - 
oiuu.u_ 920 o 


i 24+0-53 


0.56±g;« 


1220.0™ 


MWC 120 


0.868 


38^0 +140 -° 

OOOU.U_ggQ Q 


i oo+O.ll 
L ° 6 -0.16 


92+ 06 

u - yz -0.08 


1350.0™ 


VV Ser 


1.111 


2740.0™ 


1 O1+0.23 
L - J± -0.24 


o sq +0 07 

u - dy -0.07 


1200.0™ 


V1295 Aql 


1.174 


1500.0™ 


i 70+O.I2 
X - '°-0.02 


u - dz -0.01 


1270.0±g : g 


V1685 Cyg 


4.195 


171 +228ao 

1 1 1U.U_210.0 


i -- LO -o.oi 


1 -i 7 +0.05 
J--- 1 - ' -0.01 


i75o.oig:8 


AS 442 


1.187 


3580. O^ ,; , 


1 oo +0 - 12 

i - uu -0.01 


60+ 07 

u - uu -0.01 


1590.0™ 


MWC 1080 


1.224 


27060 o +5540 -° 


i yo+0.07 
x - '°-0.12 


-1 70+0. 07 
x - ' 'J-O.IS 


1740 +90 -° 



Note. - - Values of i?i nj 2, are comp uted using the deriv ed angular inner radii, 
#111,2/2, and the distances employed by Eisner et all (120041 ) . We have assumed an 
inner ring radius i?i n i=0.1 AU for these models; the high inner ring temperature 
for MWC 1080 may indicate that a larger value of i?i n ,i is warranted for this 
source ( §3.3.2j) . 
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Table 4. Best fits for disk models with temperature gradients 



Source 



T* (K) 



Xt- 



9 in /2 (mas) R in (AU) 



o 



T in (K) 



Best-Fit Disk Models 



CQ Tau 


7580 


0.976 


0.60 


f0.29 
-0.18 


09 +au4 
u - uy -0.03 


0.40 


f0.15 
-0.05 


1360.03™ 


MWC 758 


8720 


1.117 


0.63 


f0.02 
-0.12 


09 +0 ' 01 
u - uy -0.02 


0.45 


4-0.05 
-0.05 


1480.01™ 


T Ori 


8720 


0.759 


0.42 


4-0.34 
-0.32 


1Q+ - 15 
u - ±y -0.14 


0.50 


4-0.50 
-0.15 


1580.0™ 


MWC 480 


8970 


0.780 


1.03 


fO.ll 
-0.15 


14+ - 01 

u - ±4 -0.02 


0.55 


f0.10 
-0.10 


i34o.oig:8 


AB Aur 


9520 


1.657 


i 17 +0.09 
i, - l '-0.06 


16+ - 01 
u - ±o -0.01 


0.70 


f0.10 
-0.05 


1700.0^8:8 


VV Ser 


9520 


1.029 


1.00 


f0.24 
-0.22 


0-31185? 


U - DO -0.15 


1430.0^8:8 


MWC 120 


10500 


0.524 


1.06 


f0.26 
-0.22 


o.53^;i? 


0.60 


4-0.40 
-0.15 


1440.01™;° 


V1295 Aql 


10500 


1.068 


1.51 


f0.02 
-0.12 


n 44+0-01 

u -^-0.04 


1.00 


4-0.05 
-0.25 


1330.0 + ^o 


AS 442 


11900 


1.107 


0.49 


f0.26 
-0.38 


29+ - 16 

u -^ y -0.23 


0.60 


4-0.40 
-0.20 


1750.0l™ ° 


V1685 Cyg 


22000 


4.158 


0.96 


fO.Ol 
-0.04 


Q7+0- 01 

U - J ' -0.05 


1.00 


4-0.05 
-0.10 


mo.oi™ 


MWC 1080 


30000 


1.119 


0.79 


f0.10 
-0.08 


7Q+0- 11 
u - ' y -0.08 


0.65 


f0.10 
-0.05 


2060.0t?8o°o 



"Continuum" -Only Fits 



CQ Tau 


7580 


1.078 


1.16 


MWC 758 


8720 


1.218 


0.71 


T Ori 


8720 


0.648 


0.23 


MWC 480 


8970 


0.972 


0.96 


AB Aur 


9520 


0.712 


1.16 


VV Ser 


9520 


0.709 


1.03 


MWC 120 


10500 


0.737 


1.26 


V1295 Aql 


10500 


1.126 


1.54 


AS 442 


11900 


0.978 


0.10 


V1685 Cyg 


22000 


5.582 


1.00 


MWC 1080 


30000 


1.008 


0.83 



+0.04 
-0.46 
+0.18 
-0.19 
+0.51 
-0.13 
+0.24 
-0.19 
+0.19 
-0.23 
+0.21 
-0.27 
+0.04 
-0.54 
+0.03 
-0.33 
+0.34 
-0.00 
+0.01 
-0.20 
+0.20 
-0.20 



o.i7: 
o.n 



-0.01 
-0.07 
+0.03 
-0.03 

n+o- 23 

u - 11 -0.06 

n 14+0-03 

u - 14 -0.03 

16+ 003 

u - ±D -0.03 
f-0.06 
"-0.09 
H0.02 
-0.27 

4^+°- 01 
u - 4O -0.10 

06+ - 20 

u - UD -0.00 
-0.02 
-0.21 
+0.21 
-0.20 



0.32] 
0.63^ 



1.00 
0.84 



1.00 
0.50 
0.40 
0.50 
0.70 
0.70 
1.00 
1.00 
0.35 
1.00 
0.70 



0.05 
0.55 
0.15 

-0.10 
0.60 

-0.05 
0.25 

-0.10 
0.30 
0.20 
0.30 
0.20 
0.05 

-0.60 
0.05 

-0.50 
0.15 

-0.05 
0.05 

-0.35 
0.30 
0.15 



4-60.0 
-140.0 
+70.0 
-90.0 
410.0 
270.0 

i33o.otf8:8 

190.0 
140.0 
130.0 
100.0 
70.0 
210.0 
60.0 
-130.0 
80.0 
590.0 
100.0 
120.0 
220.0 
150.0 



1440.0 
1500.0 
1740.0 



1700.0 
1450.0 
1570.0 
1320.0 
2260.0 
1690.0 
2060.0 



Note. — The stellar effective temperatures, T*, are based on the assigned spectral 
types from lEisner et al.1 (J2004J). Linear inner radii, i?; n , are computed us i ng th e 
derived angular inner radii, 6- m /2, and the distances employed by lEisner et al.1 (J2004J). 
Objects are listed in order of increasing T* to illustrate the correlation of T* with i?j n , 
a, and T Ui . 
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Fig. 1. — Normalized, calibrated photometry for four Herbig Ae/Be stars and one "check 
star," HD 30111. We display these objects because the photometry for all sources was 
calibrated using the same set of three main-sequence calibrator stars, allowing examination 
of systematic calibration effects. The check star, which has a spectral type of G8III (H—K^ 
0.1 from 2MASS), exhibits a nearly flat spectrum across the if-band, illustrating that our 
calibration is reliable. The photometry for our target sources shows no evidence of narrow- 
band emission features, and the observed spectral slopes for these objects are comparable 
to those of blackbody disk models where the temperature of the hottest material is between 
1200 K (dotted line) and 1600 K (dashed line). 
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Fig. 2. — Inner ring radii, 6- m /2, as a function of wavelength, as determined from inclined 
uniform ring models fitted to the data in individual spectral channels. We als o provide the 
linea r inner ring radii, R in , computed using assumed distances for these objects (jEisner et al. 
2004 . and references therein). Inclinations, position angles, and ring widths were assumed 
( §3.2!) . and the inner ring radius was the only free parameter in the model fitting. 
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Fig. 3. — PTI V 2 data for AB Aur plotted with the best-fit single-temperature ring model 
( §3.3. ip . The model systematically under-predicts the data in the short- wavelength channels. 
The reduced \ 2 °f the fit for this object, and all other objects in our sample, is > 1, indicating 
poor agreement between this simple model and the data. 
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Fig. 4. — PTI V 2 data (points), and best-fit models. Two-ring models are represented 
by solid lines. Disk models incorporating smooth temperature gradients are indicated by 
dotted lines. Dashed lines show the best fits obtained when only the second and third 
channels ( "continuum" -only fits) were used to constrain the fits of these disk temperature- 
gradient models. Data for multiple baselines, with multiple position angles on the sky (§2]), 
are plotted here. Curves for best-fit models may therefore have a jagged appearance since 
inclined disks have different visibilities when observed at different position angles. 
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Fig. 4. — continued. 
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Fig. 4. — continued. 
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Fig. 5. — Radiative cross sections (a) for water and CO in the i^-band spectral region 
(solid line). Lines ar e drawn from the HITRAN database for temperature T = 3000 K 
(IRothman et al.ll2005l ). and convolved with Voigt profiles assuming molecular number den- 
sities of n = 10 19 cm -3 . The boundaries of the spectral channels used in this paper are 
indicated with dashed lines. 
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Fig. 6. — Best-fit values of a for disk models with a power-law temperature gradient, as a 
function of stellar effective temperature, T*. Although the uncertainties on a are large for 
some objects, the data suggest that disks around cooler stars tend to have shallower tem- 
perat ure gradients than disks around hotter stars. T* is based on the assigned spectral types 
from Eisner et al.l ((200J). For objects with the same T*, we have introduced an arbitrary 
offset of 200 K for ease of viewing. 
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Fig. 7. — SEDs computed for our best-fit models. SED data from the literature (jEisner et al. 
2004 . and references therein) are indicated with filled circles, disk models with power-law 
temperature profiles are shown with solid lines, single-temperature ring models are repre- 
sented by dotted lines, and two-ring models are shown as dashed lines. 



